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We numerically study the critical behavior of the one-dimensional XY model of the size N with 
variable interaction range L. As expected, the standard local order parameter of the magnetization 
is shown to well detect the mean-field type transition which occurs at any nonzero value oi L/N. The 
system is particularly interesting since the underlying one-dimensional structure allows us to study 
the topological excitation of the winding number across the whole system even though the system 
shares the mean-field transition with the globally-coupled system. We propose a novel nonlocal order 
parameter based on the width of the winding number distribution which exhibits a clear signature 
of the transition nature of the system. 
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I. INTRODUCTION 

The XY model is one of typical model systems that 
have been widely studied in the arena of statistical 
physics [I|. Most existing studies so far have been per¬ 
formed in the regular lattice structures of one 0 , two Q , 
three [^, and four [sj dimensions. In particular, two- 
dimensional (2D) XY model has been found to exhibit 
the famous Berezinskii-Kosterlitz-Thouless (BKT) tran¬ 
sition 3, and the all-to-all globally-coupled XY system 
has been known to exhibit the mean-field (MF) type 
phase transition Q. Also, recent progress in complex 
network research has made the interplay between connec¬ 
tion structure of interacting spins and the nature of phase 
transition a central topic in statistical physics 0,0 • Dif¬ 
ferently from most previous studies, we investigate in this 
paper the one-dimensional (ID) XY model with variable 
interaction range and explore how the interaction range 
affects the collective critical behavior of the model. The 
underlying ID topology of the system allows us to study 
the winding number excitation that has been known to 
exist in the system with local/nonlocal interaction I,®- 
We propose a novel nonlocal order parameter based on 
the distribution of the winding number excitation, which 
is found to successfully exhibit a clear signature of the 
phase transition. 

The present paper consists of five sections. Section II 
introduces the model and in Sec. Ill we explore thermo¬ 
dynamic behavior of the model, by measuring the stan¬ 
dard quantities such as magnetization, specific heat, sus¬ 
ceptibility, and Binder’s cumulant to detect the phase 
transition. In Sec. IV, which contains the key part of the 
present paper, we define the winding number for a given 
configuration of the phase variables, and investigate the 
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behavior of its distribution depending on the system size 
and temperature. Based on the observation of the wind¬ 
ing number distribution function, a novel nonlocal order 
parameter is introduced and used to successfully detect 
the phase transition. Finally, a brief summary follows in 
Sec. V. 


II. MODEL 

We study the ID XY model of N spins described by 
the Hamiltonian 

J N i+L 

i—1j—i—L 

where (j)i € (—7r,7r) is the phase angle variable of the ith 
spin with the periodic boundary condition (pi+N = <pi-, 
and L is the interaction range in one side so that each 
spin interacts with total 2L nearest-neighbor spins in 
both sides. We take the ferromagnetic coupling with 
the positive strength (J > 0), which makes the neighbor 
spins favor their phase difference to be minimized. The 
typical ID XY model with the nearest-neighbor interac¬ 
tion corresponds to the case of L = 1, and the globally- 
coupled XY model with all-to-all couplings is achieved 
when L = N/2 — 1, respectively. We in the present pa¬ 
per investigate the collective behavior of the model given 
by Eq. O, varying the interaction range L, which is a 
key control parameter that passes from the short-range 
{L/N —>■ 0) to the infinite-range (L/N —>■ 1/2) regimes 
in the thermodynamic limit of V —>■ oo. 

After suitable normalization of time, energy, and tem¬ 
perature, the first-order Langevin-type equations of mo- 
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tion of the system read 

• dH 

Vi — "TJT ^ Vi’’ 

0(pi 

^ i-\-L 

= XI sm(()!.i - (/)j)+77i, (2) 

j=i-L 

where the dimensionless thermal noise r]i satisfies 

= 2T6{t-t') (3) 

at the dimensionless temperature T in units of J/ks with 
the Boltzmann constant ks- We have taken the over¬ 
damped regime with the inertia term (containing (pi) ne¬ 
glected as in Ref. [nl where a condensed matter system 
(smectic liquid-crystal) has been studied. Equation ([2]) at 
zero temperature (T = 0) has also been used to describe 
the system of the ID coupled oscillators with variable 
interaction range If one is only interested in equi¬ 
librium behavior as in the present study, one can use the 
Mont-Carlo (MC) simulation method instead. 

We note that the ID system with short-range inter¬ 
action does not exhibit any long-range order at finite 
temperature 0. On the other hand, when the dimen¬ 
sionality of the system is higher than the lower critical 
dimension 2, the celebrated Mermin-Wagner theorem dl 
cannot disprove the existence of the long-range order, 
as typica^ exemplified in 3D Q, 4D and globally- 
coupled Q XY models. In this context, it is natural 
to have the following questions: Is there a critical inter¬ 
action range {L/N)c beyond which a hnite-temperature 
phase transition starts to occur? If so, does the phase 
transition always belong to the ME universality class? 
To answer these questions, we study the critical behav¬ 
ior of the model, first measuring some standard quantities 
for various values of L and N, which is presented in next 
section. 


III. THERMODYNAMIC BEHAVIOR 

In this section we numerically study the thermody¬ 
namic behavior of the system, with particular attention 
paid to the emergence of the phase transition at finite 
values of the interaction range L. It is possible to nu¬ 
merically investigate equilibrium behavior of the system 
in two different ways: numerical integration of equa¬ 
tions of motion Eq. (ED and MC simulation based on 
the Hamiltonian Eq. (ED- It is straightforward to show 
that both numerical methods are equivalent to each other 
as long as only equilibrium behaviors are concerned: 
the steady-state solution of the Eokker-Planck equation 
based on Eq. Q is simply the Boltzmann distribution 
P ~ exp{—H/kBT) with the Hamiltonian in Eq. ([1} [l3l| . 
Since the MC simulations usually run much faster than 
direct integration of stochastic differential equations, we 
here use the former with the standard Metropolis local 
update algorithm and measure various quantities of inter¬ 
est. In the MC simulations, all quantities are measured 
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FIG. 1: (Color online) Magnetization is plotted as a func¬ 
tion of the temperature T for various system sizes N and 
interaction ranges L for (a) L/N = 0.4 and (b) L/N = 0.1, 
respectively. 


over 10^-10® MC steps after equilibration over the initial 
10® MC steps. 

We measure the equilibrium magnetization defined as 


TO0 = / 

1 ^ 

vS'"' 

\ 

i=i 


where (• ■ ■) denotes the thermal average. The behavior of 
the magnetization is shown as a function of the tem¬ 
perature T, varying the value of L to keep the ratio L/N 
unchanged [see Eig.[T]. We find that the magnetic order¬ 
ing (m 0 > 0) occurs at T < 0.5 for L/N = 0.4, as shown 
in Fig. UDa). The signature of the transition from the 
disordered phase (m^ = 0 ) to the ordered one > 0 ) 
becomes clearer as the system size N increases. We also 
measure the magnetization for smaller value, L/N = 0.1 
[see Fig. [T](b)] and L/N = 0.05 (not shown here), and 
find that the two cases also show the magnetic order¬ 
ing for T < 0.5, for sufficiently large system sizes. The 
behavior of the magnetization shown in Fig. [1] suggests 
that the finite-temperature phase transition occurs at any 
nonzero value of L/N. In other words, we expect that the 
critical value of the interaction range beyond which the 
phase transition occurs at a nonzero critical temperature 
satisfies {L/N)c = 0 in the thermodynamic limit. 

To see it further, we also investigate other standard 
quantities such as the specific heat c«, susceptibility X; 
and Binder’s fourth-order cumulant Ub 0 , defined by 
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NkBT^ 


m - {H? 


N 


X = I (m^) - 


Ub = 1- 


(ml) 

3{mlY’ 


(5) 

( 6 ) 

(7) 


where H is the Hamiltonian in Eq. ED- Figures [5] and 
[3] show Cv and Ub as functions of T for L/N = 0.4 and 
L/N = 0.1, respectively. We find that both quantities, 
Cy and Ub, strongly support the emergence of the phase 
transition at T r; 0.5 for both L/N = 0.4 and 0.1. Al¬ 
though not shown here xC^) a peak which shifts to¬ 
wards Tc ~ 0.5 as the system size is increased both for 
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FIG. 2: (Color online) (a) Specific heat c„ and (b) Binder’s 
cumulant Ub versus the temperature T for various system 
sizes N. The interaction rage L is chosen to yield the same 
value of the ratio L/N = 0.4. 
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FIG. 3: (Color online) (a) Cv and (b) Ub versus T for L/N = 
0.1. For comparison, see Fig. [2]for L/N = 0.4. 

L/N = 0.4 and 0.1. We also investigated for the case of 
L/N = 0.05 (not shown), and obtained the same conclu¬ 
sion if we focus on larger system sizes. 

We turn our attention to the universality class of the 
phase transition, and consider the critical behavior of the 
magnetization characterized by 

( 8 ) 

with the critical exponent /3 and the critical temperature 
Tc in the thermodynamic limit. According to the finite- 
size scaling theory [Tsj . we expect that in a finite¬ 
sized system satisfies the scaling form 

m^ = N-P/-f{{T-T,)N^/-), ( 9 ) 

where is the critical exponent that describes the critical 
behavior of the correlation volume ^ in d dimen¬ 
sions: ^ \T — Tc\~'' [l^. The scaling function f{x) 

with the scaling variable x = [T — Tc)N^/^ has limiting 
behaviors: f{x) ~ as x —>■ 0 and f{x) ~ const, as 
X —>■ -boo. At criticality (T = Tc), the magnetization 
reduces to 

~ ( 10 ) 

To estimate the exponents /3 and v we plot as a 
function of the size N in the log-log scale (not shown), 
and measure its slope, which gives us the exponent fi/v- 
We also check the finite-size scaling relation directly by 
plotting versus (T — Ty)N^^'', controlling the 

values of fi/v and v [see Fig. SKa)]. We find that the 
scaling function f{x) converges to a constant for large x. 


FIG. 4: (Golor online) Finite-size scaling collapse of the 
magnetization m^\ m^N^^ " versus (T - Tc)N^/‘' for (a) 

L/N = 0.4 and (b) L/N = 0.1. With Tc = 0.5, P/v — 1/4, 
and P = 2 chosen, good quality of scaling collapse of data is 
achieved. 

and behaves as x^ for small x as expected. The numerical 
findings shown in Fig. HKa) is then summarized as 

/3/p =1/4 and P = 2, (11) 

which yields fi = 1/2. This result implies that the nature 
of the phase transition for L/N = 0.4 is the same as that 
of the MF transition @ . The crossing of the specific heat 
Cy at Tc shown in Fig.[21[a) and Fig.[31[a) also implies that 
the specific heat emonent a = 0, in accord with the MF 
universality class Q. We also examined the case for the 
smaller value of L/N = 0.1 [see Fig. |4l[b)] and L/N = 
0.05 (not shown), and obtained the same result, which 
suggests that the mean-field nature of phase transitions 
should be robust for any nonzero value of L/N. 

IV. WINDING NUMBER EXCITATION AND A 
NONLOCAL ORDER PARAMETER 

Recently, the emergence of twisted wave in the system 
of coupled oscillators with local/nonlocal interaction has 
been reported ilS- We note that the development of 
the twisted wave in dynamic models has the same phys¬ 
ical origin as the winding number excitation in the ID 
XY model. In this section, we measure the winding num¬ 
ber across the system in equilibrium for each sample and 
compute its probability distribution function by using a 
large-size ensemble of samples. We define the winding 
number q for a given phase configuration by 

1 ^ 

^ i=l 

where ‘mod’ denotes that the phase difference (fi+i — (fi 
is measured modulo 27r and the periodic boundary con¬ 
dition (f)N+i = 4>i is used. It is to be noted that the gap¬ 
less spin-wave-type excitation can occur without chang¬ 
ing the winding number. This reminds us the topolog¬ 
ical excitation of vortices in the conventional 2D XY 
model [U, Q. 

From the observation of the standard thermodynamic 
quantities in Sec. Ill, the system surely undergoes a sin¬ 
gle phase transition. Accordingly, we expect that the 
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FIG. 5: (Color online) Probability distribution function P{q) 
of the winding number q (a) at T = 0.5 for various system 
sizes N and (b) for a given size N = 640 at various tem¬ 
peratures T. For each P{q) the width of the distribution 
(7{N, L, T) is computed and then used to scale the horizontal 
axis as shown in (c) and (d), each obtained from (a) and (b), 
respectively. For (a)-(d), L/N = 0.4 has been used. 

existence of the phase transition should also alter the 
pattern of the winding number excitation in some way 
at the observed critical temperature Tc ~ 0.5. We mea¬ 
sure the probability distribution function of the winding 
numbers to see a signature of the phase transition, after 
sufficient equilibration procedure. Figure [SKa) shows the 
distribution P[q) for various system sizes N at the crit¬ 
ical temperature Tc(= 0.5). We find that the width of 
the distribution function P[q) systematically changes, as 
the system size N is increased. The inversion symmetry 
P{—q) = P{q) is easily understood since the system has 
no reason to prefer clockwise winding to counterclockwise 
one (and vice versa). The temperature dependence of the 
distribution P{q) is also investigated at a given system 
size N = 640 [see Fig. [5][b)]. Again, we observe that the 
width of the distribution changes as T is varied. In or¬ 
der to check the possibility of the scaling of the winding 
number distribution function, we measure the standard 
deviation cr defined by 

a^^{q^)-{q)\ (13) 

where {q^) = f q^P{q)dq and {q) = f qP{q)dq. We then 
scale the horizontal axis in Fig. [SJa) and (b) by using 
the scaling variable q/a with a numerically computed 
for given values of N, L, and T. It is shown very clearly 
that after the winding number q is scaled by the width a 
of the distribution function, all the curves are put on top 
of each other as displayed in Fig. [5l[c) and (d), obtained 
from Fig. [SKa) and (b), respectively. Our observation of 
the collapse of the winding number distributions strongly 
suggests that the width a{N,L,T) of the distribution 
function can successfully represent the whole distribution 
function and thus will be used below for further analysis 
to study critical behavior in the system. 


FIG. 6: (Color online) (a) Standard deviation a is plotted 
as a function of the temperature T for various system sizes 
N. (b) The newly defined nonlocal order parameter rriq in 
Eq. (O versus T for various values of N. (c) versus N in 
log-log scale clearly exhibits the power-law decay form at T = 
Tc ~ 0.5, with the slope corresponding to jdqlvq = 0.49. (d) 
Finite-size scaling collapse of mq yields Vq = 2.00(5). From 
(c) and (d), we estimate Tc = 0.500(3), Pq/vq = 0.49(4) and 
Vq = 2.00(5). All results in (a)-(d) are for L/N = 0.4. 

In Fig. [ni[a) we show cr as a function of T for various 
system sizes with L/N = 0.4 fixed. From the central- 
limit theorem for the independent random variables, we 
expect that P{q) at infinite temperature has a ~ . In 

order to compensate such a behavior, we define the nor¬ 
malized width s = a/'/N, and introduce a new quantity 
mq defined by 

mq = Sao - s, (14) 

where Soo = s{N —>• oo,T —>• oo). From Fig.[6l[a), it is ex¬ 
pected that mq —>■ 0 from above in the high-temperature 
limit and mq > 0 in the low-temperature regime, hope¬ 
fully playing the role of the order parameter. The cal¬ 
culation of the value of Soo is straightforward since the 
phase difference of the nearest neighbors (/i+i ~ 4‘i sim¬ 
ply becomes an independent random variable due to the 
lack of any spatial correlation in the high-temperature 
limit. In other words, Xi = mod {4>i+i — (f'i) inde¬ 
pendent from each other and randomly distributed in 
Xi G (—7r,7r]. From the symmetry of the distribution we 
get {q) = 0, and the second moment is computed as 

NT o N 

where P{x) is the uniform probability distribution func¬ 
tion for X and the cross terms XiXj {i ^ j) has given 
null contribution from the independence between Xi 
and Xj at infinite temperature. Consequently, we get 
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= \/( 9 ^) ~ ( 9 )^ = /(^v^) in the high-temperature 

limit, yielding Soo = o/y/N = l/(2-\/3) ~ 0.288675. We 
then plot Fig. El^b) for niq versus T, which looks very 
much similar to the standard order parameter in 
Fig. [H Finite-size scaling ansatz then allows us to ex¬ 
pect the existence of the scaling function given by 

mq = (16) 

where g(x) is a scaling function that behaves as g{x) ^ 
as a; —>■ 0, and g{x) ^ const, as a: —>■ -l-oo. At critical¬ 
ity (T = Tc), the new nonlocal order parameter niq is ex¬ 
pected to show the power-law behavior: mq ~ . 

On the basis of the prediction, we detect the exponents 
Pq and Vq^ by plotting ruq as a function of N for various 
T, where the slope at gives the value of Pqjvq [see 
Fig. IHKc)]. If only the large system sizes (iV > 160) 
are used, we find that the curves follow the power-law 
form in Fig. |ni[c) at T = 0.497-0.502 with the slopes 
Pq/^q = 0.45-0.52. For temperatures outside of this 
range, the curves exhibit clear deviations from the power- 
law form. The best fit to the power-law form is obtained 
at T = 0.500 and Pqjvq = 0.49, and thus we conclude 
Tc = 0.500(3) and Pqjvq = 0.49(4). The exponent Vq 
is obtained from the data collapse: Figure [6] (d) shows 
the behavior of mqN^‘‘/'^'> against (T — where 

Pq/vq = 0.49 and i>q = 2.00 at Tc = 0.500 are used for 
the best collapse, yielding Pq = 0.98(11). It is particu¬ 
larly important to recognize that our new nonlocal or¬ 
der parameter based on the topological winding number 
excitation gives the correlation exponent Vq = 2, identi¬ 
cal to P = 2 previously confirmed for the standard local 
magnetization order parameter vurj). However, we believe 
that further study is required to understand the value of 
Pq ~ 0.98 which is very close to unity. We also exam¬ 


ine results for L/N = 0.1 and L/N = 0.05 (not shown 
here); although not as clear as for L/N = 0.4, we are 
able to confirm the same values of exponents as long as 
sufficiently bigger system sizes are used. 


V. SUMMARY 

In summary, we have numerically investigated the crit¬ 
ical behavior of the ID XY model of N spins with vari¬ 
able interaction range L. It has been confirmed that the 
critical interaction range beyond which the phase tran¬ 
sition starts to occur at a nonzero finite temperature is 
very small, and presumably {L/N)c = 0. The nature 
of the transition has been examined by measuring stan¬ 
dard quantities such as the magnetization, specific heat, 
susceptibility, and Binder cumulant. All measured quan¬ 
tities unanimously suggest the transition is of the mean- 
field type at any nonzero value of L/N . The underlying 
one-dimensional topology of the system makes it possi¬ 
ble to study the winding number excitation. By system¬ 
atically examining the probability distribution function 
of the winding number excitation, we have suggested a 
novel nonlocal order parameter based on the width 
of the winding number distribution function. We show 
that our new order parameter rUq can successfully detect 
the phase transition and its nature. 
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